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ABSTRACT: 

A technique using hybrid Finite Element Method (FEM) /Method of 
Moments (MoM), and Geometrical Theory of Diffraction (GTD) is presented to analyze 
the radiation characteristics of cavity fed aperture antennas in a finite ground plane. The 
cavity which excites the aperture is assumed to be fed by a cylindrical transmission line. 
The electromagnetic (EM) fields inside the cavity are obtained using FEM. The EM fields 
and their normal derivatives required for FEM solution are obtained using (1) the modal 
expansion in the feed region and (2) the MoM for the radiating aperture region(assum- 
ing an infinite ground plane). The finiteness of the ground plane is taken into account 
using GTD. The input admittance of open ended circular, rectangular, and coaxial line 
radiating into free space through an infinite ground plane are computed and compared 
with earlier published results. Radiation characteristics of a coaxial cavity fed circular 
aperture in a finite rectangular ground plane are verified with experimental results. 
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1. INTRODUCTION: 


Cavity backed aperture antennas are useful in aerospace applications due 
to their nature to be conformal to the surface of the host object[l]. To predict the input 
impedance and radiation pattern of these types of antennas, one may use the Method of 
Moments. However, application of MoM to these problems result in full dense matrix 
and hence for complex cavity structure MoM will require prohibitively large computer 
storage and CPU time. Alternatively, Finite Element Method (FEM) has been more 
attractive for the analysis of these problems because FEM has a relatively simple formu- 
lation for complex penetrable structures. Also, use of FEM results in sparse, banded 
matrices, which can be effectively stored and solved. However, FEM for unbounded 
problems when used alone does not incorporate the Sommerfield radiation condition. 
The radiation conditions at the boundaries of discretized region must be imposed by 
using appropriate absorbing boundary conditions. Recent efforts have concentrated on 
the use of absorbing boundary conditions close to the radiating antenna to reduce the 
discretization region[2], Unfortunately, the accuracy of these approximate boundary 
conditions depends on specific problems, leading to results of unpredictable accuracy. 
To eliminate the disadvantages of MoM and FEM methods, a combined approach using 
FEM method inside the cavity and MoM outside the cavity has been proposed [3,4,5] by 
several authors. 

Using the combined FEM/MoM techniques, the input impedances of cav- 
ity backed microstrip antennas have been calculated in [1], The modelling of coaxial feed 
in [1] is done by replacing the feed probe by a constant filament of current with a delta 
gap generator. The feed modelling in [1] is expected to give good results for shallow 
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cavities. However, for deep cavities rigorous feed modelling is required. In this paper a 
combined Finite Element Method(FEM) and Method of Moments(MoM) approach with 
detail feed modelling is used to predict characteristics of cavity backed antennas in infi- 
nite ground plane. Various feed lines like rectangular, circular and coaxial waveguides 
exciting aperture through a cavity are considered. The electromagnetic fields inside the 
cavity are calculated using the FEM formulation. The boundary values of the field and 
their normal derivatives required in the FEM solution are obtained by using the 
waveguide modal expansion function in the feed waveguide region and in the aperture 
region by using a global method such as the method of moments. Using the present 
method input admittance of rectangular, circular, and coaxial apertures in an infinite 
ground plane and excited, respectively, by rectangular, circular, and coaxial transmission 
lines are calculated and compared with earlier results. 

In most antenna measurements, an infinite ground plane is simulated by a 
relatively large finite ground plane. The change in the radiation pattern of an aperture 
antenna due to a finite ground plane can be determined using the Geometrical Theory of 
Diffraction(GTD) [6]. In the implementation of GTD, the radiated fields from the aper- 
ture antenna are first determined assuming an infinite ground plane. These fields are 
used to calculate the incident fields on the edges of finite ground plane. Employing the 
compact edge diffraction coefficients given in [7] and the incident field on the edges, the 
diffracted field in the forward as well as backward directions is calculated using the pro- 
cedure given in [6]. The total radiation field is then obtained by superposition of far field 
with infinite ground plane assumption and the diffracted field. To verify the validity of 
the GTD technique when used in combination with FEM/MoM technique, the radiation 


3 



patterns of a coaxial cavity with finite ground plane are calculated and compared with 
experimental results. Numerical data for input return loss of such a coaxial cavity is also 
computed and compared with experimental data. 

2. THEORY 

The geometry of the problem to be analyzed is shown in figure 1. The elec- 
tric field inside the cavity satisfies the vector wave equation 

Vx^i-Vxlj-rel = 0 (1) 

where (i r , £ r are the relative permeability and relative permittivity of the medium in the 
cavity. The electric field E is solved via hybrid FEM/MoM technique assuming an infi- 
nite ground plane. Following the procedure described in [1,3]/ equation (1) can be writ- 
ten as 


III = j( 0 [i o \\[Txh]*H ap ds 

V V s a 

+ | ( T x h ) • Hinpds (2) 

s inp 

where S a is the aperture surface in the infinite ground plane and 5 ; is the common sur- 
face between input feed line and the cavity as shown in figure 1. H ap is the magnetic 
field over the aperture S a and H nif) is the magnetic field over the surface S ( . T is the 
vector testing function. 
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In order to solve equation (2) for E , the volume of the cavity is subdivided 
into small volumetric tetrahedral elements. Inside a tetrahedral element the electric field 
may be expressed in terms of the vector edge basis functions[5]. 
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Selecting the testing function T to be identical to the basis functions, the 
left hand side of the equation (2) may be written in a matrix form as 



( 3 ) 


^ is the column matrix with 

unknown coefficients for vector basis functions as its elements. 

The integeral over the radiating aperture S a is evaluated by finding the 

electromagnetic field outside the aperture by replacing the aperture with magnetic cur- 
— ^ 

rent M , backed by an infinite ground plane[l]. The discretization of cavity volume into 
tetrahedral elements automatically results in discretization of the aperture into triangu- 
lar elements at z=0 plane. The integral is then evaluated over all the triangles on the 
aperture surface using the procedure outlined by Rao et al[ 8]. It can be written in a 
matrix form as 


where the matrix z is a symmetrical sparse matrix, and 



( 4 ) 


where 


is full dense matrix and 


is a column matrix whose elements are the 


unknowns amplitudes of electric field over the aperture. 


Hinp required in the evaluation of the second term on the right hand side 
of equation (2) can be determined using waveguide vector modal functions. For a 
waveguide feed, the transverse magnetic field in the feed line can be formulated as the 
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incident field and reflected fields as 



where h m (x, y) is the waveguide vector modal function as defined in [9], a m is the 
reflection coefficient at the feed line and cavity interface, J m is the propagation constant, 
Y is the modal admittance, and H, n is the incident field. With the substitution of equa- 
tion (5), the second term on the right hand side of equation (2) can be written in matrix 
form as 

b. ] - fz. I • \e. 1 (6) 

inpj |_ inpj |_ inpj v ' 

For a detailed derivation of equation (6), the reader is referred to [10]. Using equations 
(3), (4), and (6), equation (2) may be written as 


[zl ■ Tel - [zl ■ Tel + fz. 1 • \e. 1 = \b. 1 (7) 

Lrd L5J |_ a] |_ aj |_ inpj |_ inpj |_ inp ' 

Since the column matrices e a and e jnp are the part of all Z matrices can be com- 
bined. The equation (7) then can be written as 

[a] [<?] = b inp (8) 

where ^ is a partly sparse (due to FEM) and partly dense (due to MoM and mode 
matching), is the excitation vector and f^j is the unknown vector to be solved. 

The biconjugate gradient algorithm is used with diagonal preconditioning [11] to effi- 
ciently solve equation (8). This algorithm also makes use of the symmetry of the matrix 
to effectively store the matrix. 
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Using the equation (5) the reflection coefficients of any waveguide mode 


can be obtained. 



(9) 


The input admittance is then calculated as 


Y. 


i + r 

l-r ° 


(10) 


where Y is the characteristic admittance of the feed transmission line. 

O 

Once the electric field in the aperture is known, the magnetic field at the far 
field can be evaluated directly from the equation: 


0 , <|>) 



jk sin0 (xcos(|) + y sintf)) 
(00 + <j)<j>) • M (x, y) e 


dxdy 


( 11 ) 


Effect of finite ground plane : 

The infinite ground plane solution described in the above sections yields 
the field radiated in the forward direction. This field is obtained by equation (11). The 
effect of finite ground plane is incorporated by employing the compact edge diffraction 
coefficients [7] and following the procedure described in [6]. The radiation fields are sup- 
plemented by the additional field contributions arising from the diffraction by the edges 
of the finite ground plane. The diffracted field supplements the field in the forward 
region whereas it accounts totally for the field in the shadow region(behind the finite 
ground plane) where the infinite ground plane solution gives a zero field. 

The geometry of the finite ground plane is illustrated in figure 1. The mid 
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points on the two edges of length W 2 are given by A and B, respectively, whereas the 
mid points on the two edges of length Wj are given by C and D, respectively The dif- 
fracted fields at A and B are given by[6] 


- \" S‘ U 5- °A (p -4- A PV n) 


'cA 


(r+P cA ) 


1/2 


exp (- jkr ) exp {jk ( d sin0) (12) 


h\ b ~ \H^d, f. *)d, <p cB , y, V, PV,. n) 


'cB 


r(r+p.) 


1/2 


exp (-jk r) exp {-jk d sin0)(13) 


for 4> polarization, and 


0 }D h (p c „, 5,, Vo , P'„, «) 


cA 


(r+P cA ) 


1/2 


cn/i (~jk o r) exp (Jk Q d sin0) (14) 


As ~ H/i* I’ ” J D » A ^ ") 


As 


r(r + p M ) 


1/2 


exp {-jk r) exp {-jk dsin0X15) 


for 0 polarization. In the above equations. 


and 


IT, 

P cA ~ P cB ~ d ~ ~2~ 


(16) 


t >l = ^ + 0 0 < 0 < 7t 


5-e 

o<0<5 

i 

CD 

^ < 0 < 71 
2 


(17) 


(18) 


Similarly, the diffracted fields at C and D are also calculated by replacing 

w 2 

P c c = P cD = d = ~2 • The detailed expressions for diffraction coefficients D y//; are 



given by [7, eq.113]. Total diffracted field is given by 

( — ) d — y d 

->d I Ha + Hb 
H = 1 

\ d — ^ d 

H C +H d 

The diffracted fields for (f> = 1 8(/ 
lated in a similar manner. The total radiated far field is given by adding the diffracted 
field and farfield magnetic field given in equation (25); i.e., 

—yt y -yd 

Hjf= H ff +H (20) 

It should be noted that in the above analysis, the interaction between the 
aperture field and the diffracted field is ignored. It can be seen in the following section, 
that even with this approximation, the present analysis predicts the effect of finite 
ground plane on the radiation pattern with reasonable accuracy. 

3. NUMERICAL AND EXPERIMENTAL RESULTS 

To validate the present analysis various numerical examples are consid- 
ered. As a simple case, open ended waveguide structures excited by a dominant mode 
field and radiating into half space through an infinite ground plane are analyzed. To 
study the effect of finite ground plane, a coaxial aperture in a finite ground plane and 
excited by a coaxially fed cavity is considered. The numerical data for these cases is 
compared with earlier published results and the experimental data obtained at NASA- 
Langley Research Center. 


cf) = 0 plane 
(f> = 90^ plane 


(19) 


and (f> = 270° planes can also be calcu- 
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3.1 Apertures in Infinite Ground Plane:: 

3.1.1 Rectangular Aperture: 

A rectangular aperture as shown in figure 2a with dimensions a = 4.755 cm, 
b=2.215 cm, and L= 1.9 cm in an infinite ground plane and excited by a rectangular 
waveguide loaded with a dielectric plug(e r = 3.76) is considered(figure 2(a)). The 
dielectric plug is discretized using tetrahedral elements. The air dielectric interface 
inside the waveguide is considered as input plane. The dominant TE^ 0 mode is consid- 
ered to be incident at the plane z = -L. The input admittance of this structure is calcu- 
lated as a function of frequency and plotted in figure 2(b). The numerical results agrees 
very well with the data presented by Swift[12]. 

3.1.2 Circular Aperture: 

A circular aperture in an infinite ground plane and excited by a circular 
waveguide is considered as a second example. The geometry is shown in figure 3(a). The 
circular waveguide of length L=1X from the aperture is considered as a circular cavity 
for FEM discretization. The dominant TE 11 mode is considered to be incident at the 
plane z = -L. The input admittance is computed as a function of the radius of the 
waveguide and plotted in figure 3(b) along with the numerical data given by Bailey and 
Swift[16]. Very good agreement is seen between the results computed using the present 
technique and those of [13]. 

3.1.3 Coaxial Aperture: 

As a third example, a coaxial aperture excited by a coaxial line is consid- 
ered. The geometry is shown in figure 4(a). Again the coaxial line of length L= IX from 
the aperture is considered as a coaxial cavity for FEM discretization. Only dominant 
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TEM mode is considered for both incident and reflected fields in the coaxial feed line as 


the higher order modes do not seem to be contributing much to the final result. The 
input admittance is computed as a function of frequency and plotted in figure 4(b). The 
results are compared with the available data in the literature [14]. 

3.2 Apertures in Finite Ground Plane: 

To verify the validity of the current method when combined with GTD, a 
coaxial cavity with a finite ground plane is considered. The geometry of the cavity is 
shown in figure 5. The cavity is fed by a coaxial line with the center conductor extending 
upto the radiating aperture. The cavity volume is discretized using 4474 tetrahedral ele- 
ments, resulting in 4904 unknowns. For experimental verification, the cavity with 
dimensions as given in figure 5 is fabricated and its input return loss and radiation pat- 
terns as a function of frequency are measured at NASA-Langley Research Center. The 
input return loss computed using the present method assuming infinite ground plane is 
presented in figure 6 along with the experimental results. It can be noted that the finite 
ground plane in this case does not effect the input return loss properties. 

The far field patterns using expression (30) for the structure shown in fig- 
ure 5 are calculated for two frequencies,^ = 5.0 GHz and f 2 = 9.0 GHz and are shown 
in figures 7-8 along with the experimental data taken at NASA-Langley Research Cen- 
ter's low frequency Antenna Test Facility. Good agreement between the results obtained 
by the present analysis and the experimental data supports validity of present analysis. 
Since radiation pattern of an antenna depends upon the aperture field it is informative to 
know the tangential electric field over the aperture of the structure shown in figure 5. In 
figures 9 the tangential electric field over the aperture is plotted for = 5.0 GHz and 
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/ 9 = 9.0 GHz. As expected the aperture field deviates from the TEM type field as the fre- 
quency of operation increases. 


4. CONCLUSIONS 

Radiation characteristics of cavity backed aperture antennas in finite 
ground plane are presented using a hybrid FEM/MoM technique and GTD. Various 
numerical examples are considered to demonstrate the flexibility and validity of the 
method. Very good agreement is seen between the numerical data obtained using the 
present method and the data available in the literature. The successful application of the 
GTD method in conjunction with FEM-MoM is demonstrated for calculating the radia- 
tion patterns of a coaxial cavity with finite ground plane. Since FEM is used for cavity 
region this method is applicable to arbitrarily shaped cavities with complex and aniso- 
tropic materials. By proper modelling, the present method can be used to analyze 
microstrip antenna arrays residing in a cavity with finite ground plane. This method can 
also be easily extended to a finite circular ground plane without much effort. 
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LIST OF FIGURES 


Figure 1 Geometry of a cavity backed aperture in finite ground plane. 

Figure 2 (a) Rectangular waveguide with a dielectric plug radiating into half space. 
fl=4.755 cm., £>=2.215 cm., L=1.9 cm., £^.=3.76 
(b) Variation of normalized admittance of the above structure with frequency. 
The solid and dashed lines are computed using the present analysis and 
the solid circles are from Swift[12], 

Figure 3 (a) Circular waveguide with radius r radiating into half space. L=1.0?i, E =1.0 
(b) Input admittance of the above structure as a function of 2 r/X . The solid 
and dashed lines are computed using the present analysis and the solid cir- 
cles are from Bailey and Swift[13] 

Figure 4 (a) Coaxial line with inner radius r x and outer radius r 9 radiating into half 
space. r 2 /r L =1.57,L=1.0X, £ ; .=1.0 

(b) Input admittance as a function of k r l . The solid lines are computed from 
the present analysis and the solid circles are from Swift[14] 

Figure 5 Geometry of a coaxial cavity with finite rectangular ground plane. W 1 =24", 
W2=12", L=3/8". Outer radius of the coaxial cavity=l". Inner radius of the 
coaxial cavity=0.0181". The cavity is fed by a 5011 coaxial line. 

Figure 6 Return loss of the coaxial cavity with frequency. The solid line is computed 
using the present analysis and the solid circles are the experimental results. 

Figure 7 (a) Radiation pattern of the coaxial cavity at 5 GHz at (f> =0° and (]) = 180° 

(b) Radiation pattern of the coaxial cavity at 5 GHz at (f> =90° and (f> = 270° 

This method with GTD, This 

method with infinite ground plane, _ _ Experiment. 
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Figure 8 (a) Radiation pattern of the coaxial cavity at 9 GHz at (f> =0° and (]) = 1 80° . 

(b) Radiation pattern of the coaxial cavity at 9 GHz i cf) = 90° and (]) = 270° . 

This method with GTD, This method 

with infinite ground plane, _ _ Experiment. 

Figure 9 The vector plot of the electric field at the aperture of the coaxial cavity. The 
amplitudes have been equalized for clarity of presentation. 

(a) 5 GHz (b) 9 GHz. 
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